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ABSTRACT 

In  this  paper  we  present  the  novel  of  two  immiscible  fluids  flow  and  the  effect  of  heat  and  mass  transfer  along 
a oppositely  moving  vertical  plates.  The  vertical  channel  is  subjected  to  the  transverse  magnetic  field  and  both  plates  are 
moving  with  constant  velocity  w0  in  opposite  directions.  The  coupled  governing  equations  of  the  flow  and  heat  transfer 
with  appropriate  boundary  conditions  are  solved  using  Finite  Element  method.  The  profiles  of  velocity >,  Micro  rotation 
and  temperature  are  studied  for  various  parameters  Grashof  Number  (Gr),  magnetic  field  (H),  Reynolds  Number(R), 
Eckert  Number  (Ec),  Material  Parameter  (K)  and  represented  graphically.  The  Nusselt  number  and  shear  stress  values 
are  also  analyzed. 
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INTRODUCTION 

The  subject  of  two-fluid  flow  and  heat  transfer  has  been  extensively  studied  due  to  its  importance  in 
chemical  and  nuclear  industries.  The  design  of  two-fluid  heat  transport  system  for  space  application  requires 
knowledge  of  heat  and  mass  transfer  processes  and  fluid  mechanics  under  reduced  gravity  conditions. 
Identification  of  the  two-fluid  flow  region  and  determination  of  the  pressure  drop,  void  fraction,  quality  reaction 
and  two-fluid  heat  transfer  coefficient  are  of  great  importance  for  the  design  of  two-fluid  systems.  Lohrasbi  and 
Sahai  [3]  studied  two-phase  MHD  flow  and  heat  transfer  in  a parallel  plate  channel  with  the  fluid  in  one  phase 
being  electrically  conducting.  Malashetty  and  Leela  [7]  have  analyzed  the  Hartmann  flow  characteristics  of  two 
fluids  in  horizontal  channel.  The  study  of  two-phase  flow  and  heat  transfer  in  an  inclined  channel  has  been  studied 
by  Malashetty  and  Umavathi  [5]  and  Malashetty  et  al  [6]  .Micropolar  fluids  are  non-Newtonian  fluids  with 
microstructures  such  as  polymeric  additives,  colloidal  suspensions,  liquid  crystals,  etc.  Eringen  [1]  developed  the 
theory  of  Micropolar  fluids,  in  which  the  microscopic  effects  arising  from  the  local  structure  and  the  micro 
motions  of  the  fluids  elements  are  taken  into  account.  The  study  of  viscous  dissipation  is  applicable  to  polymer 
technology  involving  the  stretching  of  plastic  sheets.  The  Problems  of  Micropolar  fluid  flow  between  two  vertical 
plates  (channel)  are  of  great  technical  interest.  A lot  of  attention  has  been  given  by  many  researchers.  Suresh  babu 
et.al  [2]  studied  the  heat  transfer  of  micropolar  and  viscous  fluids  in  a vertical  channel.  Sastry  and  Rao  [11]  have 
studied  the  effect  of  suction  in  the  laminar  flow  of  a Micropolar  fluid  in  a channel.  Bhargava  and  Rani  [8]  have 
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examined  the  convective  heat  transfer  in  Micropolar  fluid  flow  between  parallel  plates.  Stamonkovic  et.al  [9]  studied  the 
heat  transfer  of  two  immiscible  fluids  between  horizontally  moving  plates. 

Keeping  in  view  of  the  wide  area  of  practical  importance  of  multi  fluid  flows  as  mentioned,  the  objective  of  this 
study  to  investigate  the  effect  of  heat  and  mass  transfer  in  a vertical  channel  with  oppositely  moving  plates. 

Mathematical  Formulation 

The  two  infinite  parallel  plates  are  placed  at  Y=  - hi  and  Y=  h2  along  Y-direction  initially  as  shown  in  Figure  1 and 
both  plates  are  isothermal  with  different  temperaturesT , and  T2  respectively.  The  distance  from  (-hi  to  0)  represents  region 
1 and  distance  from  (0  to  h2 ) represents  region  2 where  the  first  region  is  filled  with  Micropolar  fluid  and  the  second  is 
with  viscous  fluid.  The  fluid  flow  in  the  channel  is  due  to  buoyancy  forces.  The  transport  properties  of  both  fluids  are 
assumed  to  be  constant. 


Region- 1 

Region -2 

^ 

> 

w0 

- w0 

-hi 

h2 

0 

Figure  1:  Physical  Configuration 

We  consider  the  fluids  to  be  incompressible  and  immiscible  and  the  flow  is  steady  laminar  and  fully  developed. 
The  governing  equations  are 
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To  solve  the  above  system  of  equations,  we  use  the  boundary  and  interface  conditions  proposed  by  T.Arimen  et  al 
[10]  as  follows 

Ul  - w0  at  Y - -J\  U2  — -w0  at  Y - hj 

Ui(  0)=U2(0) 

(Jul+K)^  + Kn  = jU2C^L  = 0 at  Y = 0 
V 1 ’ dY  2 dY  dY 


n-0  atY  — -hj 


For  the  corresponding  temperature  boundary  conditions  it  is  assumed  that  the  temperature  and  heat  flows  are 
continuous  at  the  interface. 

T = 7]  at  Y = ,T  = T2  at  Y = h2 


7i(o)=r2(o) 


dY  2 dY 


atY  = 0 


C = C,  at  Y = ,C  = C2  at  Y = 
C1(0)=C2(0) 


Df^L  = D ^2,  atY  = 0 
1 dY  2 dY 
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We  assume  that 
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By  introducing  the  following  non  dimensional  variables: 
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Solution  of  the  Problem 

The  coupled  governing  equations  are  solved  numerically  using  the  regular  Galerkin  Finite  Element  method  as 
given  by  J.N.  Reddy  [4].  For  computational  purpose  each  region  is  divided  into  100  linear  elements,  each  element  is  3 
nodded. 

jih 

The  shape  functions  at  each  node  of  a typical  1 element  are  the  Langrange’s  interpolation  polynomials  given  by 
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The  stiffness  matrix  equations  corresponding  to  the  governing  equations  (10)  to  (15)  for  1 element  are  evaluated 
by  using  the  following  equations: 


Region  1 


s„dui  ^du’s;  cis>’ 


dY 


+ - 


dY  dY 

1 \Gr_ 

R 


dY  + 


kl 


1 + kl 


N{S{Sp-\N{Sj-^-dY 

111  J ‘ ' dY 


l + kl 

d(N!S{)  {dSj(d(N!Sj) 


\0iS!S?dY  + 


Gc  hD-Y 


dY 


1 f 

J dY 


dY 


dY- 


R 1 + hD 

2 K1 


lS'dr  -T^lu»s'dr=0 


1 + k 


2 + Kl 


2 1 {N/S/ ) Sf’dY  - Sf  {UdS/ ) + j — ((/,/V ) dY 


= 0 


t( dS j djAjSj) 


1 dY 
Region  2 


-J 


V 


dY  dY 


dY  + Pr£cJ 


dUdSJ 


i V 


1;  i' 


Sp  dY  = 0 


5." 


jy 


-j 


jy  dY 


Y+ — 

hp 


— Uejs; ) sp  dY+-^-—Ui  -y)  s,pdy 

R J V 2‘  ' ’ 1 1 + hD  R T ’ 1 


M2m 

h 


\s'(u,js;)dY=o 


rf(02/s/) 

r 

)aS'/  d 

dY 

J 

£2 

dY 

\ 

jy 

j 

dY  + VxEc—[S; 
ml 


a rAdUJSJ 


i\ 


2 i^i 


dY 


dY  -Q 


Q & C2  are  pnear  functions  of  Y . 
Where  ^ is  the  typical  element  region 
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These  coupled  governing  equations  are  solved  iteratively  subject  to  the  boundary  conditions  given  in  (14)  until  the 
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The  results  are  depicted  in  the  graphs  from  Figure  2 to  Figure  19.  The  velocity  profiles  are  from  Figure  2 to  Figure 
7.  The  velocity  is  more  in  the  mid  region.  The  velocity  is  more  in  the  viscous  region  for  all  variations  than  in  micropolar 
region.  The  thermal  buoyancy  retards  the  flow  (Figure  2).  From  Figure  3 it  is  evident  that  the  increase  in  inertial  force 
retards  the  flow  in  micropolar  region  significantly,  but  the  variation  in  viscous  region  is  negligible.  The  Lorentz  force 
dominates  the  flow  from  Figure  4. The  heat  dissipation  is  inversely  proportional  to  the  flow.  The  effect  of  Ec  is  more 
significant  in  viscous  region  than  the  micropolar  region.  The  effect  of  material  parameter  is  minute  on  the  velocity  (Figure 
6). 


The  variation  of  Micro  rotation  (N)  with  various  parameters  is  exhibited  from  Figure  8 to  Figure  13.  The  angular 
is  momentum  is  reversal  for  all  the  variations  due  to  the  walls  movement.  The  gradual  variation  of  N is  observed  with  Gr 
(Figure  8),  R (Figure  9),  H (Figure  10),  Ec  (Figure  11),  K (Figure  12)  and  Gc  (Figure  13).  The  Lorentz  force  is  more 
predominant  in  enhancing  the  angular  momentum  where  as  all  other  forces  considered  are  inversely  proportional  to  the 
angular  momentum. 

All  the  temperature  profiles  are  exhibited  from  figure  14  - 19.  The  effect  of  Gr  (Figure  14),  Gc  (Figure  19), FI 
(Figure  16)  and  Ec  (Figure  17)  on  temperature  is  found  significant.  The  buoyancy  enhances  the  temperature.  The  Lorentz 
force  retards  the  temperature.  Obviously  the  heat  dissipation  enhances  the  temperature  and  it  is  more  pronounced  for  Ec 
>0.003.  The  effects  of  Reynolds  number  and  the  micro  rotation  parameters  on  temperature  are  not  much  significant. 


The  Nusselt  number  and  the  shear  stress  is  given  in  Table  1.  The  buoyancy  reduces  the  heat  transfer  rate  in 
micropolar  fluid  boundary  and  enhances  at  the  viscous  boundary.  The  reverse  effect  is  observed  for  shear  stress.  The 
Lorentz  force  enhances  the  heat  transfer  rate  at  the  micropolar  fluid  boundary  but  reduces  at  the  viscous  boundary.  The 
viscous  dissipation  is  also  dominant  in  heat  transfer  rate  and  shear  stress.  Material  property  reduces  the  stress  on  the 
boundary. 


Figure  2:  Variations  of  Velocity  with  Gr 


Figure  3:  Variations  of  Velocity  with  R 
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Figure  4:  Variations  of  Velocity  with  H 


Figure  6:  Variations  of  Velocity  with  K 


Figure  8:  Variations  of  Micro  Rotation  with  Gr 


Figure  10:  Variations  of  Micro  Rotation  with  H 


Figure  12:  Variations  of  Micro  Rotation  with  K 

Impact  F actor (JCC):  5.6934 


Figure  5:  Variations  of  Velocity  with  Ec 


Figure  7:  Variations  of  Velocity  with  Gc 


Figure  9:  Variations  of  Micro  Rotation  with  R 


Figure  11:  Variations  of  Micro  Rotation  with  Ec 


Figure  13:  Variations  of  Micro  Rotation  with  Gc 

NAAS  Rating:  2.45 
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Figure  14:  Variations  of  Temperature  with  Gr 


Figure  16:  Variations  of  Temperature  with  H 
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Figure  15:  Variations  of  Temperature  with  R 


Figure  17:  Variations  of  Temperature  with  Ec 


Figure  18:  Variations  of  Temperature  with  K Figure  19:  Variations  of  Temperature  with  Gc 

Table  1:  Nusselt  Number  and  Shear  Stress 


Gr 

R 

H 

Ec 

K 

St-I 

St-II 

Nu-I 

Nu-II 

100 

2 

5 

0.001 

1 

36.2402 

-28.7413 

-0.395748 

-0.663513 

200 

2 

5 

0.001 

1 

42.3718 

-31.408 

-0.361873 

-0.702001 

300 

2 

5 

0.001 

1 

48.5982 

-34.352 

-0.319367 

-0.748942 

300 

1 

5 

0.001 

1 

84.3885 

-34.5898 

-0.0843623 

-0.755735 

300 

3 

5 

0.001 

1 

36.8886 

-34.4316 

-0.347081 

-0.761117 

300 

4 

5 

0.001 

1 

31.0417 

-34.5025 

-0.351712 

-0.770067 

300 

2 

8 

0.001 

1 

32.1439 

-16.1166 

-0.458378 

-0.556303 

300 

2 

10 

0.001 

1 

26.2391 

- 11.321 

-0.476314 

-0.534776 

300 

2 

15 

0.002 

1 

17.2485 

-5.43871 

-0.476903 

-0.5424 

300 

2 

5 

0.003 

1 

50.6399 

-40.785 

0.185761 

-1.52488 

300 

2 

5 

0.001 

2 

42.738 

-34.4623 

-0.321857 

-0.75571 

300 

2 

5 

0.001 

3 

39.4375 

-34.5435 

-0.32035 

-0.761079 

CONCLUSIONS 

• The  velocity  increases  drastically  in  the  mid  region  and  more  in  viscous  region. 

• The  temperature  is  more  in  Micropolar  region  than  the  viscous  region. 
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• The  viscous  dissipation  enhances  the  heat  transfer  on  the  boundary  Y=  -hi  and  depreciates  the  heat  transfer  on  the 
boundary  at  Y=  h2.The  reverse  effect  has  been  observed  due  to  the  magnetic  field. 

• The  viscous  dissipation  depreciates  the  stress  on  both  boundaries.  The  reverse  effect  has  been  observed  due  to  the 
magnetic  field. 
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APPENDICES 

Nomenclature 


u,  ...  . 

N Ul  :Nusselt  number  at  ^ ^ 

1 : velocity  in  the  region  1 

Nu,  XI  , . v = 1 

2 :Nusselt  number  at  y 

u 

2 : velocity  in  the  region  2 

Sti : shear  stress  at  >’=“1 

jj 

0 : Average  Velocity 

St2  : shear  stress  at  ^ = 1 

^ : Temperature  of  the  plate  at  ^ ^ 

Greek  letters: 

^ 2 : Temperature  of  the  plate  at  ^ 

T 

: density  of  fluid  in  Region  1 

0 : Average  temperature. 

k 

^ 2 : density  of  fluid  in  Region  2 

1 : Thermal  conductivity  in  Region  1 

: Coefficient  of  Thermal  expansion  in  Region  1 

Impact  Factor(JCC):  5.6934 


NAAS  Rating:  2.45 
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k 

2 : Thermal  conductivity  in  Region  2 

: Coefficient  of  Thermal  expansion  in  Region  2 

g : Acceleration  due  to  gravity 
n : Micro  rotation  parameter 

^ 1 :Viscosity  of  fluid  in  Region  1 
ll 

K : Vertex  viscosity 

^2  : Viscosity  of  fluid  in  Region  2 

H 

G : electrical  conductivity 

0 : Magnetic  Field  Intensity 
Fr 

: Eckert  number 
:Groshof  number 

^:Spin  gradient 
y 

' 1 :Dynamic  viscosity  of  Micropolar  fluid 

^ e : Magnetic  Field  Permeability 
a : Thermal  conductivity  ratio 

^ : Reynolds  Number 

N : Micro  rotation  number 

^ : channel  width  ratio, 

TYl 

: viscosity  ratio 

^ : Thermal  expansion  coefficient  ratio 
0 : Magnetic  induction 

P : density  ratio 

Pf  :Prandtl  number 
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